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Abstract 

e + e _ annihilation data recorded by the JADE detector at PETRA were used to mea- 
sure the C-parameter for the first time at y/s = 35 and 44 GeV. The distributions 
were compared to a resummed QCD calculation which recently became available 
for this observable. In addition, we applied extended resummed calculations to the 
heavy and wide jet broadening variables, Bt and Bw, which now include a proper 
treatment of the quark recoil against multi-gluon emission with single-logarithmic 
accuracy. We further investigated power corrections to the mean values of the ob- 
servables mentioned above. In this study we considered all available e + e~ data 
between -v/5 = 35 and 172 GeV. 
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1 Introduction 



In a recent publication 0] we have presented a study of event shapes observables and 
determinations of a s using data of e + e~ annihilations at = 22 to 44 GeV recorded 
with the former JADE detector |J at the PETRA collider. This study provided valuable 
information which was not available before from e + e~-annihilations in the PETRA energy 
range. The results on a s , obtained in a similar manner as those from the experiments 
at LEP, demonstrated that the energy dependence of a s (Q) is in good agreement with 
the prediction of Quantum Chromodynamics (QCD). Evolved to the Z° mass scale, the 
results are in good agreement with those obtained at LEP, and are of similar precision. 
In addition, power corrections, applied to analytic QCD calculations of the mean values 
of event shape distributions, were found to qualitatively and quantitatively match the 
effects of hadronisation. Thus QCD could be tested without the need of phenomenological 
hadronisation models. 

Meanwhile the perturbative calculations for the jet broadening variables were improved 
by including a proper treatment of quark recoil 0]. Furthermore for the C-parameter a 
resummation of leading and next-to-leading logarithm terms to all orders of a s (NLLA) 
became available |J . Beside these advances in the perturbative description of event shape 
observables progress was made in the understanding of non-perturbative power corrections 
to the event shape observables and their mean values. In particular two-loop calculations 
of such corrections were performed |J which modify the one-loop result of the power 
correction to the event shapes by a factor (Milan factor). 

In this paper we complement our previous publication by a new (9 (a;? )+ NLLA de- 
termination of a s from C-parameter and update the determination from jet broadening 
at \/s = 35 and 44 GeV in the Sections |3] and [|. We also applied power corrections to 
the C-parameter distributions and re-investigated those for the mean values of the thrust, 
heavy jet mass and both jet broadening observables in Section |5|. We start in Section 
with a brief summary of the data samples used and draw conclusions from our results in 
Section ||. 

2 Data samples and Monte Carlo simulation 

We analysed data recorded with the JADE detector in 1984 to 1986 at centre-of-mass 
energies of 39.5-46.7 GeV and around 35 GeV. The JADE detector was operated from 
1979 until 1986 at the PETRA electron-positron collider at centre-of-mass energies of 
v/i = 12 to 46.7 GeV. 

A detailed description of the JADE detector can be found in |3[]. The main com- 
ponents of the detector were the central jet chamber to measure charged particle tracks 
and the lead glass calorimeter to measure energy depositions of electromagnetic showers, 
which both covered almost the whole solid angle of An. 

Multihadronic events were selected by the standard JADE selection cuts J/J. All 
charged particle tracks, assumed to be pions, with a total momentum of \p\ > 100 MeV/c 
were considered in the analysis. Energy clusters in the electromagnetic calorimeter, as- 
sumed to be photons, were considered if their energies exceeded 150 MeV after correction 
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year 


[GeV] 


data 


MC 


1984/85 


40-48 


6158 


14497 


1986 


35 


20 926 


25123 



Table 1: Number of events in data and in Monte Carlo detector simulation retained after 
application of the multihadron selection cuts described in the text. 

for energy deposited by associated tracks. 

Background from two-photon processes and r-pair events and from events with hard 
initial state photon radiation were removed by cuts on the visible energy E v i S = J2Ei, the 
total missing momentum p m i SS = | YlPi\ {Pi an d Ei are the 3- momentum and the energy 
of the tracks and clusters), the longitudinal balance relative to the e + e~ beam axis of 
momenta pbai = | Upll E v - ls \ and the polar angle of the thrust axis, 6^: 

• E vis > v/i/2 ; 

• £> miss < 0.3 • v/i ; 

• p ba i < 0.4 ; 

• | cos6> T | < 0.8 . 

Thus the backgrounds from 77 and r-pair events were reduced to less than 0.1% and 
1%, respectively H|. The numbers of events which were retained after these cuts for this 
analysis are listed in Table [T| 

Corresponding original Monte Carlo detector simulation data for 35 and 44 GeV were 



based on the QCD parton shower event generator JETSET 6.3 |T(J. These original Monte 
Carlo events at 35 GeV had the coherent branching for the parton shower while the 44 GeV 
events had non-coherent branching Q The main parameters used for event generation are 
given in Section |3.2j . Both samples included a simulation of the acceptance and resolution 
of the JADE detector. 



3 Experimental procedure 

From the data samples described in the previous section, the event shape distribution of 
the C-parameter was determined. The C-parameter is defined as flTf 

C = 3(AiA 2 + A 2 A 3 + A3A1) 

where A 7 , 7 = 1,2,3, are the eigenvalues of the momentum tensor 

Q al3 = T.iPiPi/\Pi\ 

1 The different treatment of coherence in these samples of simulated data has no visible influence on 
the results of this study. 
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3.1 Correction procedure 



Limits of the detector's acceptance and resolution and effects due to initial state photon 
radiation were corrected by applying a bin-by-bin correction procedure to the event shape 
distributions. The correction factors were defined by the ratio of the distribution calculated 
from events generated by JETSET 6.3 at hadron level over the same distribution at detector 
level. The hadron level distributions were obtained from JETSET 6.3 generator runs 
without detector simulation and without initial state radiation, using all particles with 
lifetimes r > 3 ■ 10~ 10 s. Events at detector level contained initial state photon radiation 
and a detailed simulation of the detector response, and were processed in the same way 
as the data. 

Next, the data distributions were further corrected for hadronisation effects. This was 
done by applying bin-by-bin correction factors derived from the ratio of the distribution 
at parton level over the same distribution at hadron level, which were calculated from 
JETSET generated events before and after hadronisation, respectively. The correction 
factors are typically of the order 10 to 20% growing large towards the 2-jet region. In the 
case of the C-parameter the correction factors are also large next to the 3-jet boundary 
which is at C = 0.75. The data distributions, thus corrected to the parton level, can be 
compared to analytic QCD calculations. 



3.2 Systematic uncertainties 

Systematic uncertainties of the corrected data distributions were investigated by modifying 
details of the event selection and of the correction procedure. For each variation the whole 
analysis was repeated and any deviation from the main result was considered a systematic 
error. In general, the maximum deviation from the main result for each kind of variation 
was regarded as symmetric systematic uncertainty. The main result was obtained using 
the default selection and correction procedure as described above. 

In detail, we considered either tracks or clusters only for the measurement of the event 
shape distributions. We varied the cut on cos9 T by ±0.1. The cut on p m i SS was either 
removed or tightened to p m i SS < 0.25 • a/s. Similarly, the momentum balance requirement 
was either restricted to pbai < 0.3 or dropped. We also varied the cut for the visible energy 
E v - ls by ±0.05 ■ In order to check the residual contributions from r-pair events we also 
required at least seven well-measured charged tracks. 

To study the impact of the hadronisation model of the JETSET 6.3 generator, the 
values of several significant model parameters were varied around their tuned values from 



Reference [|12] used for our main result. Different sets of correction factors to correct 
the data from hadron level to parton level were generated by varying single parameters 
of the JETSET generator. The variations were chosen to be similar to the one standard 
deviation percentage limits obtained by the OPAL Collaboration from a parameter tun- 



ing of JETSET at y/s = M^p [|T3[| . In detail, we investigated the effects due to parton 
shower, hadronisation parameters, and quark masses. The amount of gluon radiation dur- 
ing the parton shower development was modified by varying Alla by ±50 MeV around the 
tuned value 400 MeV. To vary the onset of hadronisation, we altered the parton shower 
cut-off parameter Q by ±0.5 GeV around the tuned value of 1 GeV. We used the full 
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observed variation of a s to reflect a variation of Qq between and 2 GeV. The width 
o"o = 300 MeV of the transverse momentum distribution in the hadronisation process was 
varied by ±30 MeV. The LUND symmetric fragmentation function, applied to hadronise 
events of up, down and strange quarks, was varied by changing the a parameter from 0.5 
by ±0.225 whereas the b parameter was kept fixed at 0.9. As a systematic variation we 
used the LUND 
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14] fragmentation function for charm 



instead of the Peterson et al. 
and bottom quarks. The effects due to the bottom quark mass were studied by restricting 
the model calculations which were used to determine the correction factors to up, down, 
strange, and charm quarks (udsc) only. Any deviation from our main result due to mass 
effects was treated as asymmetric error. 



4 Determination of a s 

4.1 Corrected event shape distributions 

After applying the corrections for detector and for initial state radiation effects we obtained 
the C-parameter event shape distributions at hadron level. In Tables ^| the corrected data 
values are listed with statistical errors and experimental systematic uncertainties. The 
mean values of the distributions are also given. Our measured results for the jet broadening 
event shape distributions can be found in Ref. 0. 

4.2 Determination of a s using 0{a 2 g )+NLLA calculat ions 

We determined a s by x 2 fits to event shape distributions of C and also Bt and By/ cor- 
rected to the parton level. For the sake of direct comparison to other published results we 
chose the so-called In (R) -matching scheme to merge the 0(a 2 s ) with the NLLA calcula- 
tions. The fits to the C-parameter distributions applied the resummation results obtained 
in For the fits to the jet broadening measures we used the improved calculation of 
Ref. [|J which includes a proper treatment of the quark recoil against an ensemble of soft 
gluons that is essential in the calculation of the jet broadening distributions. The renor- 
malisation scale factor, x^ = fi/y/s, was set to — 1 for the main result. Here, the value 
of \i defines the energy scale at which the theory is renormalised. 

The fit ranges for each observable were determined by choosing the largest range for 
which the hadronisation uncertainties remained below about 10 %, for which the x 2 /d.o.f. 
of the fits did not exceed the minimum by more than a factor of two, and by aiming at 
results for a s that are independent from changes of the fit range. The remaining changes 
when enlarging or reducing the fit range by one bin on either side were taken as systematic 
uncertainties. Only statistical errors were considered in the fit thus resulting in % 2 /d.o.f. 
larger than unity. The finally selected fit ranges, the a s results of the \ 2 fits and of 
the study of systematic uncertainties are tabulated in Tables ^ and pi and are shown in 
Figures || and ||. 

We also changed the renormalisation scale factor in the range of x^ = 0.5 to 2.0. 
We found variations larger than the uncertainties from the detector correction and the 
hadronisation model dependence. The dependence of the fit result for a s on x^ indicates 
the importance of higher order terms in the theory. 
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It should be pointed out that the improved perturbative calculation for the jet broad- 
ening resulted in a larger a s value from the fit to the data. The systematic uncertainties 
are not affected by the new calculation but the x 2 improved slightly 

We combined these new results with the a s results of our previous publication 0] re- 
placing the results obtained from the jet broadening observables. A single a s value was 
obtained from the individual determinations from thrust, heavy jet mass, total and wide 
jet broadening, C-parameter and differential 2-jet rate following the procedure described 



in References || [T5|, |16[ . This procedure accounts for correlations of the systematic un- 
certainties. At each energy, a weighted average of the six a s values was calculated with 
the reciprocal of the square of the respective total error used as a weight. In the case of 
asymmetric errors we took the average of the positive and negative error to determine the 
weight. For each of the systematic checks, the mean of the a s values from all considered 
observables was determined. Any deviation of this mean from the weighted average of the 
main result was taken SIS db systematic uncertainty. 

With this procedure we obtained as final results for a s 

a s (35 GeV) = 0.1448 ± 0.0010(stat.) +°;°^(syst.) 
a s (44GeV) = 0.1392 ± 0.0017(stat.) +°;°^(syst.) . 

The systematic errors at 35 and 44 GeV are the quadratic sums of the experimental 
uncertainties (±0.0017, ±0.0032), the effects due to the Monte Carlo modelling (IqIooSL 
iaoo2?) an d the contributions due to the variation of the renormalisation scale (taoo57, 
iaooss)- It should be noted that the modelling uncertainties due to quark mass effects 
contribute significantly to the total error. 



5 Mean Values of Distributions and QCD Power Cor- 
rections 

5.1 Power corrections 

The value of a s can also be assessed by the energy dependence of mean values of event 
shape distributions. Perturbative calculations exist for the mean values of thrust, heavy 
jet mass, total and wide jet broadening and C-parameter up to 0(a 2 s ). An observable T 
is given by the expression 



where the coefficients Ajr and were determined from the 0(a s ) perturbative calcu- 



lations [IT], [TS|, |TP|, ED]. The term —2Ajr accounts for the difference between the total 
cross-section used in the measurement and the Born level cross-section used in the pertur- 
bative calculation. The numerical values of these coefficients are summarised in Table EJ. 

In this study we corrected for hadronisation effects by additive power-suppressed cor- 
rections (l/-\/i) to the perturbative predictions of the mean values of the event shape 
observables. The non-perturbative effects are due to the emission of very low energetic 
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Observable JF 


A T 


Br 


a T 


(T) 
(M 2 H /s) 
(Bt) 
(B w ) 
(C) 


2.103 
2.103 
4.066 
4.066 
8.638 


44.99 
23.24 
64.24 
-9.53 
146.8 


-1 
0.5 
0.5 
0.25 

3vr/2 



Table 2: Coefficients of the perturbative prediction [|L7, 18, [19], |2(J and coefficients and 
parameters of the power corrections [fj to the mean values of the event shape observables. 
Note that the definition of a? does not include the 2j\4/it factor introduced in 0. 



gluons which can not be treated perturbatively due to the divergence of the perturbative 
expressions for a s at low scales. In the calculations of Reference |2l| and also || which 
we used in this analysis a non-perturbative parameter 



a (/ii 



Hi Jo 



dk a s (k) 



was introduced to replace the divergent portion of the perturbative expression for a s (y/s) 
below an infrared matching scale [ii- The general form of the power correction to the 
mean value of an observable T was first given in Reference [[H]]. It was the result of a 
one-loop calculation. In References || similar calculations have been performed at two- 
loops. It could be shown that the two-loop result modifies the one-loop result only by a 
factor Ai ph 1.8 which is known at the Milan factor. An additional factor 2/V is due to the 
different definitions of the non-perturbative parameters a e s, being the so-called "effective 
coupling" [p2l , and a as defined above. 

The two-loop result for the power correction assumes for thrust, heavy jet mass and 
C-parameter the form 



(J^ pow -) 



4CV 2M 



71 



71 




2tt \ Hi Po , 



where Cf = 4/3, while for the two jet broadening measures the correction is logarithmically 
enhanced by a factor \r(^/Qb)- We approximated Qb by \xi such that 



4CV 2M 



71 



71 



'J±) ■ In ( ^ 



a (m) - a s (y/s) 



27T V IX! Po 



Using this approximation an extra correction without the logarithmic enhancement 
was left out because it is small and strongly anti-correlated to the enhanced correction. 
Furthermore we found that the available data on the jet broadening are not yet sensitive 
to this extra correction. 
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The factor /3 = (11C^ — 2Nf)/Z in the two expressions stems from the QCD (3- 
function of the renormalisation group equation. It depends on the number of colours, 
Ca = 3, and number of active quark flavours Nf, for which we used Nf = 5 throughout 
the analysis. The term K = (67/18 — tt 2 /G)Ca — 5/9 ■ Nf originates from the choice of 
the MS renormalisation scheme. The remaining coefficient a? is listed in Table ^| for the 
event shapes considered. 

5.2 Determination of a s using power corrections 

We determined a s (M z o) by \ 2 fits °f the expression 

to the mean values of the five observables, thrust, heavy jet mass, C-parameter, and 
total and wide jet broadening, including the measured mean values obtained by other 



experiments at different centre-of-mass energies |23], |24], |25], |26|. For the central values of 
a s from the fits we chose a renormalisation scale factor of x^ — 1 and an infrared scale 
of m = 2 GeV. The x>/d.o.i. of all fits were between 0.8 {(M 2 H /s)) and 4.4 ((B T )). We 
estimated the systematic uncertainties by varying x^ from 0.5 to 2 and /// from 1 to 3 GeV. 

The results of the fits are shown in Figure |3| and the numeric values are tabulated in 
Table |^. It presents the values for a s and for ocq, the experimental errors and systematic 
uncertainties of the fit results. We consider these results based on power corrections as a 
test of the new theoretical prediction |J. It should be noted that the theoretically expected 
universality of «o is observed only at the level of 30%. 

Employing the procedure used in Section || to combine the individual a s values, we 
obtained 

a s (M zo ) = 0.1188 ±8:«j£ 

where the error is the experimental uncertainty (±0.0016), the renormalisation scale un- 
certainty (laooil) an d the uncertainty due to the choice of the infrared scale (^aooig), all 
combined in quadrature. This result is in good agreement with the world average value |27j 
of «7- a -(M z o) = 0.119 ± 0.006. 



6 Summary and Conclusions 

Data recorded by the JADE experiment at centre-of-mass energies around 35 and 44 GeV 
were analysed in terms of event shape distributions. 

The measured distributions were corrected for detector and initial state photon radia- 
tion effects using original Monte Carlo simulation data for 35 and 44 GeV. The simulated 
data are based on the JETSET parton shower generator version 6.3. The same event gen- 
erator was also employed to correct the data for hadronisation effects in order to determine 
the strong coupling constant a s . 

Our measurements of a s are based on the most complete theoretical calculations avail- 
able to date. For all observables theoretical calculations exist in 0(af) and in the next- 
to-leading log approximation. These two calculations were combined using the ln(i?)- 
matching scheme. We found the improved perturbative calculation of the jet broadening 
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to describe the data slightly better. With these calculations values of a s are obtained 
which are about 3% higher than previously. The a s values were also more consistent with 
those from other event shape observables. 

Combining the values of a s obtained in this analysis with those from our previous 
publication and using the new values obtained from jet broadening, the final values at the 
two centre-of-mass energies are 

a s (44 GeV) = 0.139 ±g;g^ 
a s (35GeV) = 0.145 t°;°o? , 

where the errors are formed by adding in quadrature the statistical, experimental sys- 
tematics, Monte Carlo modelling and higher order QCD uncertainties. The dominant 
contributions to the total error came from the choice of the renormalisation scale and 
from uncertainties due to quark mass effects. 

Evolving our a s measurements to y/s = M z o the results obtained at 35 and 44 GeV 
transform to 0.123 looos aric ^ 0.123 ^o'ooeL res P ec tively. The combination of these values 
gives a s (M z o) = 0.123 i°;°°f. 

The energy dependence of the mean values of the distributions can be directly compared 
with analytic QCD predictions plus power corrections for hadronisation effects involving 



an universal non-perturbative parameter «o J21], H- Our studies resulted in 



a s (M z o)= 0.119 j£gg 

which is in good agreement with our results from the (9(a^)+NLLA fits and also with the 
world average value. The universality of the non-perturbative parameter a is found only 
at a level of 30%. 
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0.079 


± 
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± 





015 


mean value 


0.3404 ± 0.0023 ± 0.0037 



Table 3: Event shape data at ^fs = 35 (left) and 44 GeV (right) for the C-parameter ob- 
servable. The values were corrected for detector and for initial state radiation effects. The 
first error denotes the statistical and the second the experimental systematic uncertainty. 
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1 V 7 . V 7 V 7 V 7 / 
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1 f\ till 1 1 \ 
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a + 0.225 


—0.0017 


/ a f\r\r\c\ 

—0.0009 
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on iv jt t r 

cr ? — 30 MeV 
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i r\ r\r\ 1 a 
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i on tv jt t r 
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—0.0014 


LUND symmetric 


i a nni 7 
+U.UU1 1 


l A A A 1 £ 
+ U.UU10 


l A AAOFC 
+ U.UUZO 


Q + 500 MeV 
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±0.0008 


±0.0009 


±0.0011 


MC modelling 


+0.0072 
-0.0032 


+0.0054 
-0.0026 


+0.0073 
-0.0044 


= 0.5 


-0.0100 


-0.0061 


-0.0107 


x„ = 2.0 


+0.0127 


+0.0081 


+0.0133 


Total error 


+0.0151 
-0.0112 


+0.0111 
-0.0085 


+0.0155 
-0.0120 



Table 4: Values of a s (44 GeV) derived using the 0(a^)+NLLA QCD calculations with 
— 1 and the ln(i?)-matching scheme, fit ranges and xVd.o.f. values for each of the 
three event shape observables. In addition, the statistical and systematic uncertainties are 
given. Where a signed value is quoted, this indicates the direction in which a s (44 GeV) 
changed with respect to the standard analysis. The scale uncertainty and quark mass 
effects are treated as asymmetric uncertainties of a s . 
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±0.0009 
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tracks only 


-0.0010 


-0.0006 


-0.0011 


clusters only 


-0.0009 


-0.0018 


-0.0007 


COS#t 


±0.0001 


±0.0002 


±0.0002 


Pmiss 


±0.0001 


±0.0003 


±0.0002 


Pbal 


±0.0002 


±0.0006 


±0.0002 


N ch 


+0.0005 


+0.0006 


+0.0007 


E ■ 
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LUND symmetric 
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A — 50 MeV 


—0.0021 
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A + 50 MeV 


+0.0018 
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+0.0020 


udsc only 


+0.0086 


+0.0077 


+0.0078 
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±0.0007 


±0.0008 


±0.0008 


MC modelling 


+0.0099 
-0.0048 


+0.0089 
-0.0045 


+0.0097 
-0.0058 


x M = 0.5 


-0.0107 


-0.0075 


-0.0110 


x M = 2.0 
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+0.0096 
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Total error 


+0.0169 
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+0.0134 
-0.0091 
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Table 5: Values of a s (35 GeV) derived as in Table (| but at 35 GeV. 
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(a) 


(1-T) 


(Ml/s) 


(Br) 


(B w ) 


(C) 


average 


a s {M z o) 


0.1197 


0.1140 


0.1211 


0.1217 


0.1201 


0.1188 


Q range [GeV] 


13-172 


14-172 


35-172 


35-172 


35-172 




rVci.o.f. 


43.2/24 


12.1/16 


38.9/9 


21.6/9 


11.0/7 




experimental 


±0.0013 


±0.0010 


±0.0020 


±0.0020 


±0.0016 


±0.0016 


= 0.5 


-0.0050 


-0.0026 


-0.0039 


-0.0002 


-0.0046 


-0.0023 


x M = 2.0 


+0.0061 


+0.0037 


+0.0049 


+0.0011 


+0.0057 


+0.0033 


to = 1 GeV 


+0.0026 


+0.0013 


+0.0044 


+0.0027 


+0.0025 


+0.0024 


to = 3 GeV 


-0.0020 


-0.0011 


-0.0036 


-0.0021 


-0.0020 


-0.0019 


Total error 


+0.0068 
-0.0055 


+0.0040 
-0.0030 


+0.0069 
-0.0057 


+0.0035 
-0.0029 


+0.0064 
-0.0053 


+0.0044 
-0.0034 



(b) 


(1-T) 


{Ml/s) 


(Br) 


(B w ) 


(C) 


a 


0.510 


0.616 


0.396 


0.325 


0.443 
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±0.012 


±0.018 


±0.019 


±0.020 


±0.011 


Xfj, = 0.5 


+0.004 


+0.012 


+0.005 


+0.082 


+0.006 


= 2.0 


-0.002 


-0.005 


-0.001 


-0.036 


-0.003 


Total error 


+0.013 
-0.012 


+0.022 
-0.019 


+0.020 
-0.019 


+0.084 
-0.041 


+0.013 
-0.011 



Table 6: Values of a s (M z o) (a) and do (b) derived using to — 2 GeV and x^ = 1 and the 
0(a 2 s ) calculations and two-loop power corrections which include the Milan factor ||. Fit 
ranges and % 2 /d.o.f. values for each of the five event shape observables are specified. In 
addition, the statistical and systematic uncertainties are given. Where a signed value is 
quoted, this indicates the direction in which a s (M z o) and o>q changed with respect to the 
standard analysis. The renormalisation and infrared scale uncertainties are treated as an 
asymmetric uncertainty on a s (Mzo). These uncertainties are treated equally for «o but 
exclude the infrared scale uncertainty. 
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Figure 1: Distributions measured at y/s = 44 GeV and corrected to parton level are shown 
for the total and wide jet broadening B T and B w and for the C-parameter. The fits, using 
the improved C(af)+NLLA QCD prediction for the jet broadening are overlayed and the 
fit ranges are indicated by the arrows. The error bars represent statistical errors only. 
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Figure 3: Energy dependence of the mean values of thrust (1 — T), heavy jet mass 
(Mff/s), total (B T ) and wide jet broadening (B w ), and of the C-parameter (C) are 
shown 25, 26]. The solid curve is the result of the fit using perturbative calculations 

plus two-loop power corrections which include the Milan factor || while the dashed line 
is the perturbative prediction for the same value of a s (M z o). 
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